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Statement  of  Work 

Random  variables  anting  in  communications,  radar,  photodetection,  and  sig- 
nal  processing  often  have  moment  generating  functions,  characteristic  functions, 
or  probability  generating  functions  for  which  formal  mathematical  expressions  are 
readily  derived  theoretically.  Determining  the  probability  distributions  of 
these  variables  from  generating  functions  analytically  is  seldom  possible.  The 
cumulative  distribution  of  such  a  variable,  important  because  it  provides  false- 
alarm,  detection,  and  error  probabilities,  can  be  expresssed  as  a  contour 
integral  in  the  complex  plane  whose  integrand  involves  the  generating  function. 

It  has  been  found  that  the  distribution  can  be  accurately  computed  numerically  if 
the  contour  is  taken  through  a  saddlepoint  of  the  entire  integrand.  We  propose  a 
thorough  study  of  the  potentialities  of  this  method  of  saddlepoint  integration 
for  distributions  of  both  continuous  and  discrete  random  variables  of  Importance 
in  the  above-mentioned  fields.  When  the  random  variable  results  from  a  quadratic 
functional  of  a  Gaussian  random  process,  the  generating  function  involves  the 
Fredholm  determinant  of  the  autocovariance  function  of  the  process.  Its  calcula¬ 
tion  for  points  on  the  contour  of  integrate  u  by  integrating  vector  Riccati  equa¬ 
tions  will  be  Investigated  for  processes  that  can  be  modeled  as  the  output  of  a 
linear  system  driven  by  white  noise.  A  coherent  signal  component  of  the  Input 
process  can  also  be  handled  by  this  analysis.  Acceleration  of  the  convergence  of 
the  numerical  integration  by  optimum  choice  of  the  contour  will  be  examined. 

Particular  applications  treated  in  the  proposal  are  to  the  detectability  of 
fading  radar  echoes,  the  distribution  of  the.  average  power  of  a  Gaussian  random 
process,  the  performance  of  the  optimum  and  threshold  detectors  of  a  narrow-band 
Gaussian  stochastic  signal  in  white  noise,  the  distribution  of  the.  number  of 
photoelectrons  counted  when  a  mixture  of  coherent  and  incoherent  light  falls  on 
an  ealsslve  surface,  and  the  distribution  of  the  number  of  output  electrons  from 
an  avalanche  photodiode.  Still  other  applications  will  be  considered  if  time  and 
resources  permit. 
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1.  Computing  Cumulative  Distributions  from  Moments 
This  research  is  centered  on  the  problem  of  computing  the  distribution  of  a 


random  variable  V  when  one  knows  its  moments 


bn  -  E(Vn),  n  *  0,  1,  2,  ...  . 


(1.1) 


From  these  one  forms  the  moment-generating  function  (m.g.f.) 


b(s)  =  vin(-s)n/n!  -  1  +  a^(-s)n 


■/' 


p(V)e”8V  dV, 


(1.2) 


which  is  simply  the  Laplace  transform  of  the  probability  density  function 
(p.d.f.)  p(V)  of  the  random  variable  V.  It  is  assumed  that  h(s)  is  regular  at 
least  in  a  vertical  "convergence  strip"  parallel  to  and  containing  the  Im  s-axis 
of  the  complex  s-plane. 

The  inverse  transform  to  (1.2)  is 


h(s)esV  ds/2iti. 


(1.3) 


where  c  lies  in  the  convergence  strip.  Our  concern  is  primarily  with  determining 
the  cumulative  distribution  function,  which  is  useful,  for  instance,  in 
specifying  false-alarm  and  detection  probabilities  in  radar  and  error 
probabilities  in  communication  systems.  Integrating  (1.3),  one  writes  this  as 


q"(v) 


/V  /•c+i“ 

p(V')dV'  -  /  s  h(s)e8^  ds/2ni. 


c  >  0.  (1.4) 
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The  complementary  cumulative  distribution  is 


q+(V> 


p(V')dV' 


h(s)esV  ds/2iri. 


c  >  0.  (1.5) 


These  integrals  can  seldom  be  worked  out  in  closed  form,  and  one  must  resort  to 
numerical  methods. 

The  most  common  procedure  is  to  utilize  the  Edgeworth  series.  One  defines 
the  cumulant  expansion  by 

00 

*n  h(s)  -  ^  bn(“s)n  =  -Vs  +  j  o2 s2  +  r(s), 
n*l 


r(s)  =2^  Kn(-s)n/n!, 


(1.6) 


n-3 


where  V  »  E(V),  oz  »  Var  V,  and  kq  is  the  n-th  cumulant  of  V.  When  the  m.g.f. 


h(s)  is  not  given  in  closed  form,  the  coefficients  b  =  x  / n!  of  this  series  can 

n  n 


be  obtained  from  the  coefficients  a  -  u  /n!  of  (1.2)  by  the  recurrent  relations 

n  n 


^n+1  "  an+l 


as 

-srS  ‘Va+l-r-  n>1- 


(1.7) 


r»l 


starting  with  b^  •  a^  »  V. 
One  then  expands 


h(8>  ■  exp(-Vs  o2s2 


>(i+L 


[r(s)] 


k 


k! 


(1.8) 


k-1 


writes  out  the  terms  in  the  brackets  as  a  series  in  powers  of  (-s),  substitutes 
into  (1.3),  and  integrates  term  by  term  to  obtain  a  series  for  the  density 


function  p(V).  This  series  Is  in  turn  Integrated  term  by  term  to  produce,  after 
appropriate  rearrangement,  the  well-known  Edgeworth  series.  It  Is  an 
“asymptotic"  series  In  the  sense  that  although  its  terms  at  first  decrease  in 
magnitude,  the  series  diverges  when  carried  beyond  a  certain  point.  Those  terms, 
proportional  to  derivatives  of  the  error  function,  are  oscillatory  functions  of 
their  argument  (V-V)/o,  and  the  terms  of  high  order  attain  large  positive  and 
negative  values.  Computation  with  the  series  is  difficult  to  manage,  therefore, 
and  the  results  are  unreliable  in  the  tails  of  the  distribution  where 
|V  -  V|  »  o. 

In  our  approach  one  avoids  the  expansion  in  (1.8),  which  is  the  cause  of  all 
the  trouble,  and  simply  integrates  (1.4)  or  (1.5)  numerically.  In  order  for  the 
numerical  integration  to  require  as  few  steps  as  possible,  one  would  like  to 
integrate  along  the  path  of  steepest  descent  of  the  integrand.*  This  path 
crosses  the  Re  s-axis  at  a  saddlepoint  Sq  of  the  integrand,  which  is  a  root  of 
the  equation 

♦  '(s)  h(s)  +v  -1=  °,  (1.9) 

where 

*(s)  =  fcn  h(s)  +  sV  -  In  (is)  (1.10) 

is  the  "phase"  of  the  integrand  in  (1.4)  or  (1.5).  [(±s)  =  +s  for  q”(V)  in 

(1.4);  (±s)  =*  -s  for  q+(V)  in  (1.5).]  Along  the  path  of  steepest  descent,  which 
may  have  several  branches,  Im  t(s)  is  constant.  On  the  principal  branch  crossing 
the  Re  s-axls  at  s  =  8^,  Im  4>(s)  *>  0.  In  many  problems  this  is  the  only  branch. 

Because  computing  the  path  of  steepest  descent  would  be  time-consuming,  one 
approximates  it.  The  simplest  approximation  is  a  straight  vertical  line  through 
the  saddlepoint  Sq;  that  is,  one  takes  c  ■  Sq  in  (1.4)  or  (1.5).  The  number  of 


steps  of  numerical  integration  can  often  be  significantly  reduced,  however,  by 
choosing  instead  a  parabolic  contour,  on  which 

x  -  sQ  +  j  tcy  ,  s  =  x  +  iy, 

with  the  curvature  <  approximating  that  of  the  path  of  steepest  descent. ^  The 
integrals  to  be  evaluated  can  then  be  put  into  the  form 

[h(s)eSV(+s)_1(l  -  ixy)]dy/jr, 

0 

8  ”  sQ  +  y  <y2  +  iy,  (1.11) 

with  sq  the  appropriate  root  of  (1.9).  For  V  »  V  =  E(V),  it  is  most  accurate  to 
compute  q+(V),  and  for  that  integration  8q<  0;  for  V  «  V  one  computes  q  (V), 
for  which  Sq  >  0.  For  V  near  V  it  does  not  matter  which  one  evaluates. 

One  solves  (1.9)  for  the  saddlepoint  most  expeditiously  by  Newton's  method, 
replacing  a  trial  value  of  Sq  at  each  stage  by 


q^V) 


8o  +  8o  -  TO^y* 


(1.12) 


Because  of  the  concavity  of  4>(s)  for  s  ■  x  on  the  real  axis,  this  method 
converges  rapidly.  If  it  is  inconvenient  to  calculate  the  second  derivative 
$"(s),  the  secant  method  can  be  used  instead  to  solve  $'(s)  -  0,  s  **  Sq.  If  even 
the  first  derivative  is  inconvenient  to  calculate,  one  can  solve  the  equation 

Im  $ ( Sq  +  ie)  *  0  (1.13) 

by  the  secant  method  for  a  sufficiently  small  positive  value  of  c. 


The  curvature  of  the  path  of  steepest  descent  at  the  saddlepoint  is  given 


When  the  variable  V  is  close  to  the  mean  V,  however,  the  path  of  steepest  descent 
may  begin  at  s  *  Sq  to  curve  in  one  direction  and  then  bend  back,  and  go  off  to 
infinity  in  the  opposite  direction  before  the  value  of  the  integrand  has  become 
sufficiently  small  so  that  the  numerical  integration  can  be  stopped.  This 
behavior  can  be  avoided  by  taking  the  curvature  of  the  parabolic  contour  of 
integration  Instead  as 

K  *  i  *',,(80)/*,,<so)»  u.io 

where 

*(s)  =  In  h(s)  +  Vs  (1.15) 

is  the  modified  phase,  which  dominates  the  variation  of  the  integrand  far  from 
the  saddlepoint.  Because  the  curvature  k  is  not  needed  to  great  accuracy,  it 
often  suffices  to  take  it  simply  as 

k  ■  Kj/3o^,  (1.16) 

where  ■  E(V-V)^  is  the  third  cumulant  of  V  and  ■  Var  V. 

It  sometimes  happens  that  one  knows  the  moments  vn  of  the  random  variable  V, 
but  has  no  analytical  form  for  the  ra.g.f.  h(s).  One  can  then  express  the  phase 
♦  (s)  as 

00 

•(.)  =  ^  ^  bn(-s)n  +  Vs  -  in  s,  (1.17) 

n=l 

whose  coefficients  bn  **  K^/n!  can  be  determined  from  the  quantities  an  ■  y^/nl  by 
the  recurrence  in  (1.7).  One  stops  the  summation  in  (1.17)  when  the  term  bn(-s)n 


becomes  sufficiently  small.  Alternatively  one  can  take 

$(s)  =*  £n  h(s)  +  Vs  -  £n  s 

and  compute  h(s)  from  the  power  series  in  (1.2).  In  many  problems  the  magnitude 
of  the  integrand  of  (1.4),  (1.5),  or  (1.11)  decreases  to  a  level  at  which  the 
integration  can  be  terminated  while  the  point  s  *  x  +  iy  is  still  within  the 
domain  of  convergence  of  the  series  in  (1.2)  and  (1.17). 


2.  The  Integrated  Output  of  a  Linear  Rectifier 


In  the  receiver  of  a  radar  searching  for  a  target  at  a  particular  range,  the 
output  of  a  narrowband  filter  matched  to  the  echo  signal  is  rectified  and  sampled 
at  an  appropriate  time  to  provide  a  datun  v.  The  data  v^  ,  v^,  ...»  acquired 
during  M  successive  interpulse  intervals  are  summed,  or  "integrated,"  to  yield  a 
statistic 


V  "  v!  +  v2  +  +  V  (2a) 

and  if  this  exceeds  a  certain  decision  level  Vq,  the  radar  decides  that  a  target 
is  present. 

When  the  noise  is  Gaussian  and  the  rectifier  has  a  quadratic  characteristic, 
the  detection  probability  q+(V^ )  =  Pr  (V  >  Vq)  is  given  by  the  M-th-order 
Q-function.^  We  have  previously  treated  the  computation  of  this  function  by 
saddlepolnt  integration.^  If  the  rectifier  is  linear,  on  the  other  hand,  the 
distribution  of  V  is  not  known  in  analytical  form,  Marcum^  used  the  Edgeworth 
series  to  approximate  the  false-alarm  and  detection  probabilities  of  a  system 
integrating  the  outputs  of  a  linear  rectifier.  This  problem  is  conveniently 
handled  by  the  method  described  in  Section  1. 

With  appropriate  normalization  we  can  write 


l(sx  +  xk)2  +  (s  +yk)2]1/2 


(2.2) 


where  x^  and  y^  are  independent  Gaussian  random  variables  with  mean  zero  and  unit 
12  2 

variance,  and  S  ■  (s  +  s  )  is  the  input  signal- to-noise  ratio,  assumed  the 
z  x  y 

same  for  all  k,  1  <  k  <  M.  Then  the  probability  density  function  of  the  v^'s  is 
p(v)  **  v  exp(-  jv2  -  S)  Iq(/2Sv)  ,  v  >  0,  (2.3) 


where  Iq(* )  is  the  modified  Bessel  function  (ref.  3,  p.  168).  The  moments  of  vk 


were  given  by  Rice  (ref.  6,  eq.  (3.10-12),  p.  107)  as 


Pn  =  E(vkn)  **  2n/2r(|+  1)  1F1(-n/2,  lj  -S)  (2.4) 

in  terms  of  the  confluent  hypergeometric  function 


r-1  '•r'; 


(2.5) 


When  S  *  0  the  m.g.f.  of  vk  is  given  in  closed  form  by-' 


hQ(s)  -  Etexp(-vks)]  =  1  -  (2n)1/2  s  es  /2  erfc  s. 


(2.6) 


where 


erfc  8  «= 


is  the  error-function  integral.  The  m.g.f.  of  V  is  then 

h(s)  =  [h0(s)]M. 


(2.7) 


(2.8) 


When  n  is  large,  the  n-th  moment  un  is  asymptotically 
Un  5  2n/2r(£+  l)e"S/2  I0(/2(n+l)S) 

by  eq.  (13.5.13),  p.  508  of  Abramowitz  and  Stegun.^  Using  Sterling's  formula  for 
the  factorial  and  the  asymptotic  formula  for  the  modified  Bessel  function,  we 
find  for  n  »  1 

aR  «  p^/n!  ~  (8nnS)  ^2  n-11^2  exp[— ~  +  (2nS)^2]. 

These  coefficients  decrease  so  rapidly  with  increasing  n  that  the  power  series  in 
(1.2)  for  the  m.g.f.  of  the  vk  converges  everywhere  in  the  complex  s-plane. 


The  first  four  coefficients  of  this  series 


ys)  -  y-s)k 

k=0 


(2.9) 


are 

aQ  -  1,  aj  -  (tt/2)  1/2  i,  1;-S) 

a2  -  1  +  S,  a3  =*  (ir/8)1/2  ^(-3/2, 1;-S) .  (2.10) 


The  rest  can  be  evaluated  by  the  recurrence 

2(m+14S)(m-l)a  -  ma  , 
m  ta—z 


a 


n+2 


(2.11) 


(m  -  l)(tn  +  2) 

so  that  the  series  in  (2.5)  needs  to  be  summed  only  twice.  One  can  then  use  the 
recurrence  in  (1.7)  to  determine  the  coefficients  of  the  series 


in  hQ(s)  -  y-s>k*  (2.12) 

k=0 

Multiplying  these  by  M  yields  the  terms  in  the  expression  (1.17)  for  the  phase 
♦(s)  of  the  integrands  of  (1.4),  (1.5),  and  (1.11). 

The  false-alarm  probability  Qq  «•  Pr  (V  >  Vq|  S  «  0)  can  be  computed  by 
numerical  quadrature  of  (1.11),  using  (2.6)-(2.8)  to  calculate  h(s)  on  the  path 
-  integration.  The  saddlepoint  sq  was  calculated  by  Newton's  method  as  in 
(1.12),  and  the  curvature  <  of  the  parabolic  path  of  integration  was  determined 
by  (1.14).  By  combining  this  method  with  the  secant  method  we  were  able  to 
determine  the  decision  level  Vq  to  attain  a  pre-assigned  false-alarm  probability 
Qq.  The  detection  probability  =  Pr  (V  >  |  S  >  0)  could  then  he  calculated 

by  (1.11),  in  which  the  ro.g.f.  h(s)  was  computed  as  in  (2.9)-(2. 12),  with 
(2.8).  This  method  worked  very  efficiently. 


3.  The  Average  Power  of  a  Random  Process 
Methods  for  calculating  the  cumulative  distribution  of  the  average  power 


dt 


of  a  Gaussian  random  process  z( t)  with  a  given  mean  s(t)  and  a  given  auto- 
covar lance  function  4>(t-s)  were  described  In  Sec.  3  of  our  1982-83  annual 
report.®  Research  Into  these  Is  continuing  with  emphasis  on  our  adaptation  of 
the  Kreln-Levlnson  method  in  part  (b)  of  that  section.  The  accuracy  of  the 
method  has  been  improved  by  computing  the  estimator  kernel 


Au(x,r)  =  Au(T»r)  +  u(l+ru) 

from  an  integral  equation  for  Au(t,r)  Instead  of  directly  from  eq.  (3.22)  of 

Q 

that  report,  following  the  suggestion  made  on  p.  20  thereof.  The  new  Integral 
equation  has  been  solved  numerically  by  the  method  of  conjugate  gradients  In 
combination  with  the  Sobolev  Identity  (eq.  (3.25)  of  ref.  8).  The  solution 
enables  more  precise  computation  of  the  moment-generating  function  h(u)  of  V  at 
points  u  far  out  on  the  contour  of  integration  in  (1.11),  and  the  results  for  the 
cumulative  distribution  of  V  are  consequently  more  accurate. 


4.  Posi tive- Integer-Valued  Random  Variables 
Let  V  be  a  random  variable  whose  value  can  be  only  a  positive  Integer  or 


zero.  If  pk  =  Pr  (V  =  k),  k  »  0,  1,  2,  the  probability  generating  function 

(p.g.f.)  of  V  is  defined  by 

00 

h(z)  =  E(z^)  g  ^  ^  pnz°.  (4*1) 

n=0 


In  terms  of  it  the  probabilities  pk  are  given  by  the  contour  integral 


/■ 


h(z)  z'(n+1)  dz/2ri. 


(4.2) 


where  C  is  a  closed  contour  enclosing  the  origin,  but  no  singularities  of  h(z), 
all  of  which  lie  outside  the.  unit  circle.  The  cumulative  distribution  of  V  is 
given  by 


k-1  r 

pr  (v  <  k)  -  P„  -  / 

n=0  *C 


z  ^h(z)  dz 
1  -  z  2iri 


and  its  complement  by 


(4.3) 


qk+  =  Pr  (V  >  k)  =  ^  Pn  =  f  ffr*  (4*4) 

n=k  JC+ 

Here  C_  and  C+  are  contours  like  C,  except  that  C+  encloses  the  point  z  =  1  and 
C_  does  not. 

As  described  in  Section  4  of  our  1982-83  annual  report,8  these  integrals  can 
often  be  expeditiously  evaluated  by  deforming  C_  and  C+  into  contours  passing 
through  saddlcpoints  x^  and  x^  of  their  integrands  and  approximating  the  paths 
of  steepest  descent.  The  saddlepoints  are  roots  of  the  equation 


where 


*'(z> 


h'(z)  _  n_  _  1 

h(z)  z  z-1 


(4.5) 


f(z)  =  in  h(z)  -  n  In  z  -  fcn[±(z-l)] 


(4.6) 


is  the  "phase"  of  the  integrand.  Applications  to  the  distribution  of  the  number 
of  electrons  at  the  output  of  a  photomultiplier  were  treated  in  that  report. 
Further  applications  will  now  be  outlined. 
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5.  Photoelectron  Counting 

When  light  falls  on  a  photoelec trlcally  emissive  surface,  the  probability 
generating  function  of  the  number  of  electrons  ejected  during  an  interval  (0,  T) 


is  given  by 

h(z)  =  E[e1(2_1)], 


(5.1) 


where  I  is  a  random  variable  proportional  to  the  total  energy  falling  on  the 
surface  during  that  interval.  If  we  assume  the  light  to  be  plane-polarized, 
quasimonochromatic,  and  spatially  coherent  over  the  entire  emissive  surface,  its 
field  can  be  written 

v(t)  =  Re  [V(t)e1<c],  0  <  t  <  T, 


in  which  the  complex  envelope  V(t)  is  a  circular  complex  Gaussian  random  process, 
and  fl  =  2xv  is  the  central  angular  frequency  of  the  light.  Then 


I 


dt 


(5.2) 


with  n  a  constant  directly  proportional  to  the  quantum  efficiency  of  the  emissive 
surface  and  inversely  proportional  to  the  photon  energy  hv. 

The  expected  value  E[V(t)]  >*  S(t)  represents  the  complex  envelope  of  a 
coherent  light  signal,  such  as  the  output  of  a  pulsed  laser.  The  difference 
vb(t)  *•  V(t)  -  S(t)  represents  incoherent  light  such  as  is  produced  by  an 
ordinary  incandescent  source.  It  is  assumed  to  be  a  stationary  random  process 
with  a  temporal  coherence  function  defined  by 

C(trt2)  -  Y  E[Vb(tl)V*(t2)l.  (5.3) 

The  spectral  density  of  the  light  is 


The  mean  number  .1q  of  photoelectrons  ejected  by  this  incoherent  component  is 


nn  -  nG(0)T; 


(5.5) 


the  mean  number  n_  ejected  by  the  coherent  component  is 


ng  =  y  n  /  |S( t> | 2  dt. 


(5.6) 


Q 

The  p.g.f.  h(z)  can  be  expressed  as7 


h(z)  *  lD(C)]_1  exp{ng(z--l)  /  /  o*(t)i|;(t,u;&;T)a(u)  dtdu}, 

0  *0 


K  -  nQ(l-z), 


(5.7) 


where 


D(5)  -  (l+*kO 
k=*l 


(5.8) 


is  the  Fredholm  determinant  associated  with  the  integral  equation 


P-1 


/ 


Xkfk(t)  -  T  /  Y(t-s)  fk(s)  ds,  y(t)  -  G(t)/G(0),  (5.9) 

J0 


whose  eigenvalues  are  Furthermore,  o(t)  is  a  normalized  signal, 

rT 


•If 


0(t)-s(t)//  |S(u)r  du, 

J0 

and  + ( t , u ; ^ ;T )  is  the  resolvent  kernel  defined  by  the  integral  equation 


(5. 1C) 


\|»(t,u;C;T)  +  ♦(t,v;Cjt)Y(v-u) 


CT~ly(t-u),  0  <  (t,u)  <  x  <  T. 


(5.11) 


In  the  past  the  distribution  of  the  number  of  photoelectrons  has  been 
evaluated  by  decomposing  the  field  V(t)  into  temporal  modes  by  a  Karhunen-Loeve 
expansion,  whereupon  the  p.g.f.  can  be  written  as  a  product  of  factors  involving 
the  eigenvalues  X^  of  (5.9)  and  the  coefficients 


f*(t)  S(t)  dt 


of  the  coherent  component.  The  distribution  corresponding  to  each  factor  or  mode 
has  the  Laguerre  form,  and  the  distribution  of  the  total  number  V  of 
photoelectron  counts  is  the  convolution  of  as  many  of  these  as  necessary.*® 

Other  modal  decompositions  can  also  be  used.**  (Further  references  are  to  be 
found  in  ref.  9.)  Such  methods  require  much  computation  time  and  storage, 
particularly  when  the  product  WT  of  the  counting  time  T  and  the  bandwidth  W  of 
r(u>)  is  large.  The  method  of  saddlepoint  integration  promises  much  greater 
efficiency  in  calculating  the  counting  distribution. 

When  the  light  lias  a  Lorentz  spectral  density. 


r(«)  -  ^oif 

<u  -Hi 


Y(x)  =  e 


the  Fredholm  determinant  has  the  well-known  form1 


DU)  =  e-ra[(g+l)2emg  -  (g-l)2e"in8J/(4g), 
g  «  (l+2e/m)1/2,  K  «  nQ(l-z),  ra  -  pT. 


(5.12) 


(5.13) 


If  the  signal  Is  a  sinusoid  of  constant  amplitude  and  with  a  frequency  displaced 


•"*.**  .**  .  »  .**  V*  *»'  *»*  *.*  .  *'  •  .  .V  ,•  .*  .•  .■  ■  •  »  ,  .  , 


g 

from  tJ  by  A,  as  we  have  shown,  the  p.g.f.  h(z)  is 

2 

h(z)  =  [DCS)]-1  exp{n  (z-l)[-~r- 

g +6 

,  H{(K+l(R-62)em8  -  (g-l)(a^2)e~°’g  ~  2r[(1-62)cos  m6  -  26  sin  m6]}i. 

(g+l)2em8  -  (g-l)2^8 

(5.14) 

with  6  “  A/u  and  the  rest  of  the  notation  as  in  (5.13).  After  putting  this  into 

the  integrals  (4.3)  and  (4.4)  we  have  computed  the  complementary  cumulative 

distribution  q^+  by  evaluating  the  integrals  numerically  along  a  parabolic 

contour  passing  through  the  appropriate  saddlepoint.  Even  for  values  of  k  of  the 

order  of  thousands,  accuracy  of  seven  significant  figures  was  attained  by  fewer 

than  twenty  steps  of  numerical  Integration. 

For  more  complicated  spectra  r(oi)  the  p.g.f.  h(z)  cannot  be  written  down  in 

closed  form.  It  can  be  obtained  numerically  at  points  on  the  contour  of 

g 

integration,  however,  from  the  alternative  expression 


-1 


h(z)  =  [DU)] 


x  exp{ng(z-l)  /  {o*(t;5*)-o*(t)J [o(t; £)-o(T)]dr}  (5.15) 

•'o 


where 


f- 

'(t;5)  =  /  i|<(T,u;Oo(u)du 

•'rv 


(5.16) 


and  <Kt;u;5)  *  ,Kt,u;5;t)  is  the  solution  of  the  Integral  equation 


♦  <t,u;0  +  KT 


I  «Kt,v; 

•'r* 


C)Y(v-u)dv 


(5.17) 


-  5T-iY(f-u),  0  <  u  <  T  <  T. 
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In  terms  of  this  function  the  Fredholm  determinant  is  given  by 


D(C)  -  exp[ 


I  ?(t,r;5)dT]. 

/a 


(5.18) 


The  function  4>(t,u;£)  is  the  kernel  of  the  minimura-mean-square-error  causal 
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estimator  of  the  process  (£/T)  V^*)  in  the  presence  of  circular-complex  white 

Gaussian  noise  with  unit  spectral  density.  As  such  it  can  be  computed  by  solving 
Kalman-Bucy  or  Chandrasekhar  differential  equations  whenever  one  can  establish  a 
state-space  model  for  generating  V^(t)  from  white  noise.  Indeed,  the  procedures 
for  determining  the  distribution  of  the  average  power  of  a  Gaussian  random 

O 

process,  described  in  Sec.  3  of  our  1982-83  annual  report,  can  be  applied  to 
evaluating  the  terms  in  h(z)  in  (5.15).  In  particular,  one  can  solve  the 
integral  equation  (5.17)  numerically  by  sampling  in  x  and  u,  replacing  the 
integral  by  a  quadrature  formula,  and  solving  the  resulting  linear  simultaneous 

O 

equations  by  the  method  of  conjugate  gradients  as  outlined  there. 

When  the  product  WT  of  the  counting  time  T  and  the  bandwidth  W  of  the 

incoherent  light  is  large,  one  can  utilize  the  so-called  "Toeplitz  approximation" 

q 

to  write  the  p.g.f.  as7 


h(z)  s  ex 


p{n  (z-1)  f  |21-t  f  inll+ST-1?(m)]dM/2Tt}, 

-L  1-KT  ^(m)  **  J 


r(u)  -  r(u>)/G(0),  E(ui)  -  /  0<t)e"iutdt,  e  -  n0(l-z).  (5.19) 


This  can  be  substituted  into  (A. 3)  or  (4. A)  for  numerical  integration  to 


V. 


•.*  •/  ^ 


6.  Distribution  of  the  Output  of  an  Avalanche  Photodiode 
In  an  avalanche  photodiode  primary  electrons  generated  at  one  face  by 
incident  light  are  accelerated  under  an  applied  voltage  and  create  hole-electron 


electrons  are  counted,  must  then  be  computed  by  summing  the  probabilities  pn. 


•  ■  -  ■  ► 
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(6.5) 


and  the  complementary  cumulative  distribution  is  =  1  -  .  Computing  these 

probabilities  in  this  manner  is  tedious  because  of  the  complexity  of  McIntyre's 
formula^  for  Pr  (n|  N),  and  when,  as  usually  in  a  multiplicative  device  such  as 
this,  the  numbers  k  are  large,  many  summations  are  required  and  round-off  error 
is  serious. 

The  method  of  saddlepoint  integration  has  been  found  most  efficient  for 
determining  the  cumulative  distribution  q^  of  the  number  of  output  elec¬ 
trons.^  One  computes  q^~  or  qfc+  by  numerical  integration  of  (4.3)  or  (4.4),  in 
which  now 


h(z)  -  expfX(M-l)], 


(6.6) 


M  and  z  being  related  by  (6.2).  Rather  than  solving  (6.2)  for  M  at  each  point  z 
on  the  contour  of  integration,  we  change  variables  in  (4.3)  and  (4.4)  and  use  M 
as  the  variable  of  integration,  obtaining  the  formula 


*  =  ±  f  e^(M_1 )  M-n{1  + 

* 


1  -  a  +  a(l-b)M  i  dM 


M  -  {I  +  a(M-l)] 


1  ari 

b'J  2iri  ’ 


(6.7) 


— 

where  C  and  C+  are  curves  In  the  M-plane  passing  through  the  image  Mq  or  M^  of 
the  saddlepoint  x^-  or  Xq+  of  the  integrand  in  (4.3)  or  (4.4).  The  M-plane  is 
cut  along  the  real  M-axis  from  a  point 


M.  -  (l-a)/[a(b-l )]  >  1 


to  M  *  ",  corresponding  to  a  cut  in  the  z-plane  along  the  Re  z-axis  from  a  point 
x^  to  »,  in  order  to  render  the  solution  of  (6.2)  for  M  single-valued.  A 
parabolic  contour  in  the  cut  M-plane  was  found  to  yield  to  seven  significant 
figures  with  fewer  than  twenty  steps  of  numerical  integration,  even  for  values  of 
k  as  high  as  50000.  Indeed,  the  larger  the  gain  G  and  the  mean  number  X  of 
primary  electrons,  the  more  accurate  this  method  of  saddlepolnt  integration 
becomes . 

In  an  optical  receiver  utilizing  an  avalanche  photodiode,  Gaussian  noise  is 

added  to  the  diode  output  in  a  subsequent  amplifier.  The  performance  of  such  a 

receiver  was  calculated  by  Personick,  Balaban,  Bobsin,  and  Kumar***  by  convolving 

the  distribution  of  the  number  of  output  electrons,  as  computed  by  McIntyre's 

formula,*^  with  a  Gaussian  distribution  representing  the  noise.  This  computation 

was  so  lengthy  that  they  resorted  Instead  to  Monte  Carlo  simulation  and  a  crude 

Gaussian  approximation.  This  problem  should  yield  easily  to  the  method  of 

saddlepolnt  integration.  The  moment-generating  function  (m.g.f.)  of  the  sum  V  of 

the  number  of  output  electrons  and  Gaussian  noise  with  mean  zero  and  variance 
2 

o*  is 

h(8>  «  E[e  sV]  =  exp[X(M-l)  + o2 s2],  (6.8) 

with  8  “  e”z  and  z  related  to  M  again  through  (6.2).  The  cumulative  distribution 
of  the  noisy  output  V  and  its  complement  can  be  computed  by  numerical  integration 
of  (1./+)  and  (1.5).  Once  again  we  expect  to  achieve  this  by  using  M  instead  of  s 
as  the  variable  of  integration. 

Extending  a  previous  analysis,*2  we  have  applied  the  saddlepolnt 
approximation10  to  this  problem,  evaluating  the  probabilities  q  (V)  and  q  (V)  as 

q*(V)  s  (2ii't"(s0)]  exp  $(s0),  (6.9) 
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where  the  phase  $(s)  is,  as  in  (1.10) 


*(s)  =  X(M-l)  +y  o2s2  +  Vs  -  in  (is),  (6,10) 

and  the  saddlepoint  Sq  is  the  root  of 

*'(s)  aX^-+o2s  +  V-  -=  0  (6.11) 

ds  s 

lying  to  the  left  of  the  origin  for  q+(V)  and  to  the  right  of  the  origin  for 
q  (V).  By  using  M  as  the  independent  variable,  (6.11)  could  be  solved  by 
Newton's  method. 

In  the  binary  receiver  analyzed  by  Personick  et  al.,^  the  decision  level  Vq 

above  which  it  decides  that  a  signal  is  present  is  set  to  minimize  the  average 

probability  of  error.  As  0's  and  l's  were  taken  as  equally  likely,  this  is  the 

point  at  which  Pq(V)  =  Pj(V),  where  Pq(V)  and  pj(V)  are  the  probability  density 

functions  of  V  under  the  hypotheses  Hq  (no  signal)  and  (signal  present), 

respectively.  These  density  functions  can  also  be  calculated  approximately  by 

1  ft 

the  saddlepoint  method  ; 

P0(V)  s  hQ(s0)  exp(s0V)[2Ti$|J(s0)]"li/2, 

Pl(V)  5  h^sj)  exp(siV)[2Ti$’i'(s1)l"1/2,  (6.12) 

where  h^(s)  is  the  m.g.f.  of  V  under  hypothesis  and 

4>^(s)  =  in  h^(s)  +  sV,  i  =  0,  1, 

is  the  associated  modified  phase.  The  saddlepoints  s^  and  Sj  are  the  roots  of 
the  equations 


4> '  (s)  =  h.  (e)  +  V  *■  0 


i  -  0,  1, 


(6.13) 


which  can  be  solved  by  a  few  iterations  of  Newton's  method.  One  then  solves  the 
equation 


In  Pl(VQ)  -  in  p0(VQ)  =  0 

for  Vq  by  the  secant  method,  starting  at  an  initial  trial  value  halfway  between 
the  means  GX^  and  GX^  of  the  output  V  under  the  two  hypotheses.  Personick  et 
al .*6  determined  the  mean  number  Xj  of  primary  electrons  needed  to  attain  an 
average  error  probability  of  10-^.  ^or  this  one  applies  the  secant  method  taking 
that  mean  number  as  the  independent  variable.  By  using  the  saddlepoint 
approximation  in  (6.9)-(6.13)  we  were  able  to  reproduce  the  results  plotted  In 
ref.  16  in  about  half  an  hour's  computation  on  an  HP-71  portable  computer. 


7.  Significance  Probabilities  for  Bank  Tests 
The  method  of  saddlepoint  integration  described  in  Sec.  4  has  been  applied 


which  equals  the  number  of  times  an  x  precedes  a  y  in  a  random  arrangement  of  the 

M  x's  and  N  y's  according  to  their  values.  Under  the  null  hypothesis  they  are 

all  independent  and  identically  distributed.  These  statistics  form  the  basis  of 

frequently  used  statistical  hypothesis  tests  known  as  "rank  tests." 

Table  of  the  cumulative  distributions  of  W  and  V  or  U  and  tables  of  the 

percentage  points  of  these  distributions,  used  in  assessing  the  significance  of 

the  outcome  of  such  a  rank  test,  take  up  a  great  many  pages,  even  for  values  of  M 
20 

and  N  up  to  50.  Although  recurrent,  relations  exist  whereby  the  distributions 
can  be  computed,  the  algorithms  are  extremely  time-consuming  and  require  a  large 
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amount  of  computer  memory  when  M  and  N  are  large.  One  would  like  to  have  a  means 
of  computing  the  cumulative  probability  or  its  complement  for  the  integer¬ 
valued  random  variables  W  and  V  (or  U)  directly. 

The  p.g.f.'s  of  these  statistics  under  the  null  hypothesis  are^® 


Vz)  =  E(zW)  =  Y\ 


(7. A) 


Vz)  =  E(zU)  “  c  (— -j-— ),  c  *=  (N^)-1.  (7.5) 

•  «  z  —1 


Because  of  the  symmetry  of  the  distributions  of  W  and  U  about  their  means, 


E(W)  =  M(M+l)/4,  E(U)  =  MN/2, 


(7.6) 


it  suffices  to  compute  the  probabilities  Pr  (W  <  k)  and  Pr  (U  <  k)  for  integers  k 
up  to  the  mean  value  of  the  statistic,  and  this  can  be  done  by  numerical 
integration  of 


f  V>  u 

J  l-z  2* 


Pr  (W  <  k) 


(7.7) 


with  a  similar  formula  for  Pr  (U  <  k).  For  the  latter,  no  generality  is  lost  by 
taking  M  <  N.  The  integrand  has  a  single  saddlepoint  Xq  on  the  Re  z-axis  in 
0  <  Re  z  <  1.  The  contour  C  was  taken  as  a  vertical  straight  line  through 
z  -  xq,  completed  by  the  portion  of  the  unit  circle  lying  to  the  left  of  it. 
Extensive  numerical  comparison  with  the  exact  probabilities  showed  that  for 
M  £  30  the  contribution  of  the  circular  part  of  the  contour  C  to  the  integral  is 
negligible,  and  it  suffices  to  approximate  the  probability  in  (7.7)  by 


*  t.  •• 


2  “  x0  +  iy 


(7.8) 


Pr  (W  <  k)  s 


The  integrand  of  (7.8)  drops  off  rapidly  as  y  increases  from  0,  and  for  M  of  the 
order  of  50  or  more  it  is  unnecessary  to  carry  the  integration  all  the  way  to  the 
unit  circle  (y  =  Y).  The  probability  Pr  (U  <  k)  is  computed  in  the  same  manner. 

The  amount  of  memory  required  by  this  method  —  in  contrast  to  algorithms 

previously  used  —  is  independent  of  M  and  N.  It  was  even  possible  to  program  it 

on  an  HP-15C  calculator.  The  amount  of  computation  time  is  roughly  proportional 

to  M  because  of  the  need  to  compute  the  p.g.f.  in  (7.4)  or  (7.5)  at  each  point  on 

the  contour  of  integration.  This  saddlepoint-integration  method  enables  the 
significance  probability  of  a  set  of  ranked  data  to  be  computed  directly  and 
eliminates  the  need  for  cumbersome  tables  of  distributions  and  percentage  points. 


8.  Robust  Radar  Detectors  (J.  A.  Ritcey) 

Mean-level  detectors  (MLD)  are  commonly  used  In  radar  to  maintain  a  constant 
false-alarm  rate  (CFAR)  when  the  background  noise  level  Is  unknown.  In  our 
annual  report  for  1982-83  we  reported  on  the  application  of  saddlepoint 
integration  techniques  to  radar  detection  problems  where  the  test  statistic  is 
compared  to  either  a  fixed  threshold  or  to  an  adaptive  threshold  composed  of  a 
weighted  estimate  of  the  mean  level  of  the  background  noise,  derived  from  a  given 
reference  channel.^  The  method  was  developed  for  an  arbitrary  number  of  pulses 
noncoherently  integrated  and  for  both  a  nonfluctuating  and  a  chi-square  fluctu¬ 
ating  target.  Recently,  we  have  considered  a  variant  of  the  basic  MLD  that  is 
more  robust  in  the  presence  of  interference  in  the  reference  channel,  and  we  have 
obtained  exact  analytical  results  for  the  probabilities  of  detection  and  false 
alarm  when  the  interference-to-noise  ratio  (INR)  is  nonzero.  In  practice,  the 
interference  could  result  either  from  returns  from  real  objects  such  as  other 
targets  or  clutter,  or  from  pulsed-noise  jamming.  Here  we  consider  only  Swerling 
II  target  returns  and  only  single-pulse  processing  (no  noncoherent  integra¬ 
tion).  The  approach  is  to  write  the  detection  probability  as  a  contour  integral, 
which  is  evaluated  by  the  method  of  residues.  Variants  of  this  system  in  which 
noncoherent  integration  is  employed  or  other  target  models  are  considered  would 
require  numerical  contour  integrations.  In  the  model  we  consider  here, 
developing  an  analytical  expression  for  the  moment-generating  function  (m.g.f.) 
of  the  reference-channel  estimate  is  the  primary  mathematical  problem. 

Radar  detection  of  a  known  signal  in  additive  white  Gaussian  noise  of 
unknown  variance  is  often  accomplished  by  comparing  the  test  statistic  for  a 
single  range  cell  with  an  adaptive  threshold  equal  to  a  scaled  estimate  of  the 
unknown  noise  variance.  For  a  system  that  quadra tically  rectifies  the  output  of 


a  matched  filter  to  obtain  the  test  statistic,  the  problem  can  be  modeled  by  the 
following  hypothesis  testing  problem: 


V  V  -  i  (x2  +  y2) 

Hi:  v  “  T  ((x+®)2  +  <y+t>2);  Ml2  ■  s2  +  t2*  (8.1) 


The  random  variables  x,  y  are  Independent  Gaussian  noise  samples  with  zero  mean 
and  unit  variance.  The  target  echo  provides  a  fixed  signal- to-noise  ratio  (SNR) 
for  a  steady  target,  and  the  amplitude  can  be  taken  as  a  random  variable  when  the 
echo  fluctuates.  The  phase  of  the  complex  amplitude  d  *  s  +  it  is  uniformly 
distributed  over  (0,  2ir),  and  the  total  SNR  is 


S  -  |dj2/2. 


(8.2) 


Implementing  the  generalized  likelihood  ratio  test,  the  system  decides  for  Hq  or 
Hj  according  to  whether 


v 


> 

< 


ar 


H 


0 


(8.3) 


where  r  is  an  estimate  of  the  noise  variance  and  a  is  a  positive  threshold  param¬ 
eter.  Typically  the  radar  uses  the  n  range  cells  surrounding  the  cell  under  test 
to  compute  an  estimate 


r 


(8.4) 


where 


1  ,  2  .  2  v 

2  (xJ  +  *J> 


(8.5) 


and 


the  Xj*  yj  J 


1,2,  ...,  n  are  independent  and  identically  distributed 


(l.l.d.)  zero-mean  Gaussian  random  variables,  and  we  have  Incorporated  a  scaling 
factor  n“*  Into  the  threshold  parameter  a.  In  a  locally  homogeneous  noise 
environment  the  x^,  y j  all  have  common  unit  variance.  The  performance  is 
completely  specified  by  the  false  alarm  and  detection  probabilities 


Qo 


(v)  dv  dr 


(v)  dv  dr 


(8.6) 


(8.7) 


where  f(r)  is  the  probability  density  function  (p.d.f.)  of  the  estimate  r  and 
p1(v),  i  *  0,  1,  is  the  p.d.f.  of  v  under  each  hypothesis.  An  equivalent 
hypothesis  test  Is 


2  =  v  - 


(8.8) 


where  the  random  variable  z  has  a  p.d.f.  pz(z)  whose  Laplace  transform  Is  given 
by 


hz(u) 


f 


-U2  , 

e  P2(z)dz 


'  f  e~uv  p(v)dv  f 

/ft  /ft 


earu  f(r)dr 


(8.9) 


or 

h  (u)  «*  <e  ux>  =  <e  uvXeuar>  =  h(u)  d(-au)  (8.10) 

z 

where  h(u),  d(u)  are  the  moment-generating  functions  (m.g.f.)  of  v  and  r 
respectively.  If  we  specify  a  Swerling  II  target  model,  the  m.g.f.  is  given  by 


h(u)  *•  (1+bu)-*;  b  ■  1  +  S 


(8.11) 
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and  S  =  y  <|d|  >  is  the  average  SNR.  The  estimate  r,  in  a  homogeneous 


environment,  has  the  m.g.f.  d(u)  given  by 


d(u)  =  (1+u) 


(8.12) 


The  detection  performance  can  be  written  as  contour  integrals  involving  these 


m.g.f. 's.  Since  h_(u)  and  p_(z)  are  related  by  an  inverse  Laplace  transform,  we 


have  directly  from  (8.8)  that 


Qd  ®  Pr  (z  >  0 


H'’  “// 


euz  h(u)  d(-au)dz 


(8.13) 


Q  *  -  /  u"1  h(u)  d(-au) 

I  2iti 


(8.14) 


where  the  contour  of  integration  C_  consists  of  a  vertical  path  in  the  complex 


u-plane  crossing  the  negative  real  axis  at  u  *  c^»  b  <  c^  <  0.  The  contour  is 


closed  in  an  infinite  semi-circle  in  the  left  half-plane  (LHP).  Since  d(u)  is 


the  m.g.f.  of  a  positive  random  variable,  the  only  singularities  lie  in  the 


LHP.  Therefore,  d(-au)  is  analytic  in  the  LHP  and  (8.14)  can  be  immediately 


evaluated  in  terms  of  the  residue  at  the  simple  pole,  u  =  -b“*.  Thus,  the 


detection  probability  is  simply  related  to  the  m.g.f.  of  the  estimate  r  by 


Qd  -  d(a/b);  b  -  1  +  S, 


(8.15) 


while 


Qq  “  d(a), 


The  key  point  of  this  result  is  that  given  a  closed-form  analytical  expression 


for  the  m.g.f.  d(u)  we  immediately  have  a  result  for  the  detection  performance. 


u  .... 


s  Vi 
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Therefore,  we  are  encouraged  to  investigate  other  systems  that  develop  a 

different  estimate  r,  but  whose  m.g.f.  can  still  be  obtained. 

To  reduce  the  degradation  caused  by  interfering  targets  falling  within  the 

*>  1 

reference  channel,  Rickard  and  Dillard  proposed  a  class  of  detectors  2>A, 

1  -  0,  1,  2,  ...,  n-1,  where  the  l  largest  range  cells  are  censored  from  the 
noise-variance  estimate.  If  the  quadra tically  rectified  output  for  the  j-th 
range  cell  is  denoted  by  Zj,  the  detector 3^  uses  the  unsealed  estimate 


(8.16) 


where  k  «  n-1  and  0  <  z^  <  <  ...  <  z^  <  «  are  the  order  statistics  of 

the  sample  (z^,  z^ . zr).  Since  r  is  a  biased  estimator  of  the  unknown  noise 

variance  we  consider  the  more  general  statistic 


r  -Z)  i(j)  +c2(k) 


(8.17) 


for  some  constant  c  >  0.  The  contour  integral  (8.14)  requires  the  m.g.f.  of  r, 
d(u). 

For  a  homogeneous  noise  environment  where  the  Zj's  are  l.i.d.  with  p.d.f. 
p(z)  *■  e“z,  z  >  0,  the  m.g.f.  of  r  can  be  shown  to  be  given  by 


J-l 


(8.18) 


From  (8.18)  we  note  that  the  noise  level  estimate  will  be  unbiased  if  we  scale  r 
by  k“*  atld  choose  c  ■  n  +  1  -  k.  Then 


k  "  k  Cl)  z(1)  +  (u+1~k)z(k)l 


(8.19) 


Is  an  unbiased  estimate  of  the  unknown  noise  variance  that  is  robust  with  respect 
to  the  presence  of  several  large  interference  pulses  in  the  reference  samples 

(zl»  **•»  zn^* 

The  performance  of  the  CMLD  in  a  multi pie- target  environment  can  be  obtained 
directly  from  (8.15),  given  the  m.g.f.  d(u).  It  can  be  shown  that  the  joint 
m.g.f . 


hf(u,s)  =  <exp[-u(z^j  +  ...  +  ^)-szm]>, 


(k-1)'  ~(k)J 


(8.20) 


where  z,,.  <  <  ...  <  z/( ,  <  ...  <  z,  .  are  the  order  statistics  drawn  from 

(1)  (2)  (k)  (n) 

an  independent  and  exponentially  distributed  sample  (z^,  Zj,  ..*,  z^  ...»  zn)  and 
where 


<z1> 


a  *  1=1,  ...,  I 

1  i  =  l*fl ,  . » . ,  n 


cf 1  =  1  +  INR 


(8.21) 


is  given  by 


where 


and 


i-1 
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(a-1 )(i+r ,(i))  +  0.-1 
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EL-  E  Ei* 
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t^l  *2=1 


Vi=1  ** 


Ll  *  Z2*  *  l*-i 


♦j  =  n  +1  -  j 

=  s  +  u  (k  -  j) 
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The  multiple  sun  in  the  second  line  of  (8.22)  is  present  only  when  i  <  l.  The 
reason  we  have  evaluated  the  joint  m.g.f.  h, (u, s)  is  twofold:  First,  the  CMLD 
with  arbitrary  bias  constant  c  has  a  reference  channel  m.g.f.  given  by  h^(u,cu) 
when  i  interferers  are  present.  Second,  we  can  evaluate  the  performance  of  the 
so-called  Order  Statistic  Detector  (OSD)  described  by  Rohling  ,  which  was  the 
estimate  r  *  2(k)»  a  single  order  statistic.  This  system  may  be  considerably 
easier  to  implement  than  the  CMLD. 

As  an  example  of  a  particular  result.  Figure  8.1  shows  the  additional  SNR 
required  (often  called  the  CFAR  loss)  for  the  MLD,  CMLD,  and  OSD  in  a  multiple 
target  environment  to  achieve  a  given  at  fixed  Qq  over  that  required  by  the 
Neyman-Pearson  detector,  which  is  the  optimal  detector  when  the  background  noise 
level  is  known.  Notice  that  although  the  CFAR  loss  increases  with  1NR  for  the 
MLD,  the  CMLD  and  OSD  have  a  bounded  loss  as  the  INR  ♦  ».  These  results  will  be 
more  fully  reported  in  a  paper  in  preparation  for  the  IEEE-AES. 
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